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Abstract 

We prove that the Black Saturns are stably causal on the closure 
of the domain of outer communications. 



1 Introduction 

The Black Saturn solution of Elvang and Figueras is a five dimensional black 
hole with disjoint event horizons with topology IxS 3 and K x S 1 x S 2 [HE]- 
It describes a spherical black holaS surrounded by a black ring. 

In this article, we address the question of causality violations in the Black 
Saturn solution. The absence of closed causal curves is one of the desired 
properties of a solution to the Einstein equations. Such property should 
be stable against small perturbations. Therefore, the closed causal curves 
are not allowed in any Lorentzian metric that is sufficiently near the original 
one. This leads to the notion of stable causality as introduced by Hawking [3]. 
Technically, a spacetime (M, g^) is stably causal if and only if there exists 
a differentiable function / on M such that V^/ is a past directed timelike 
vector field [2E]- The non- unique function / may be interpreted as a cosmic 
time that increases along every causal curve. 

The natural guess for the Black Saturn is that the generalized Weyl t 
coordinate is a cosmic time in the domain of outer communications (d.o.c). 
Under this hypothesis, the problem of stable causality of the d.o.c. was re- 
duced in [1] to the study of the component of the metric. Namely, stable 
causality will result from the following inequality in the generalized Weyl 
coordinates 

g(Vt,Vt) = g tt = - 9 -^ < , (1.1) 



1 The adjective "spherical" refers to the topology of the horizon. 
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where G y is a non-negative function which is zero only on the axis of the 
rotation of the Killing field d^. In other words, it is sufficient to show that 
vanishes as fast as G y on the axis of the rotation of and that > 
in the remaining part of the d.o.c. 

Numerical evidence for positivity of g^ in the plane of the ring (as defined 
in [21 Section 3.8]) was already signalled in the original paper of Elvang and 
Figueras |2J. This numerical evidence was extended in [T] to the part of 
spacetime covered by the generalized Weyl coordinates away from points 
where vanishes. Before this work, all analytical proofs of stable causality 
were restricted to the situation when the Komar angular momentum of the 
spherical component of the horizon is equal to zero: under this restriction, 
the positivity of g^ in the plane of the ring was established in [2] , and stable 
causality of the d.o.c. was shown in [1]. 

In this article, we prove that the d.o.c. of the general Black Saturn solution 
is stably causal. In addition, we show that the event horizons may be included 
to the domain of stable causality, hence the Black Saturns are stably causal 
on the closure of the d.o.cJl 

The calculations presented her^] involve manipulations of huge algebraic 
expressions and were done with Mathematica. We use the same notation 
and definitions of auxiliary functions as in pp. 

2 Stable causality 

In the generalized Weyl coordinates (t,ip,<p, p, z) the d.o.c. corresponds to 

{p > 0} U {p = 0, z [a 5 , a 4 ] U [a 3 , a 2 }} , 

where a^i,...^ are parameters. Since g^ = g^(p, z,a±, . . . ,a§), then we 
would like to prove stable causality for a union of the following sets in IR 7 

2 In fact, our proof remains valid for a family of the Black Saturn solutions with conical 
singularities on the axes of the periodic Killing fields d^, d v . 

3 The Mathematica code is available at |http : //th . if . uj . edu . pl/-szy bka/BScausality 
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V a = {ai < a 5 < a 4 < a 3 < a 2 } , (2.1) 

Vj = { P >o}nv a , 

V n = {p = 0,z<a 1 }nV a , 

Vin = {p = 0,a 1 <z<a B }nV a , 

Viv = {p = 0, a 4 < z < a 3 } D V a , 

% = {P = 0, a 2 < z}HV a . 

Hence, we have Vd. . c . = V U V}/ U V/// U V/y U V^. We are interested in 
non-degenerate solutions, so the parameters are restricted to V a . The "plane 
of the ring" corresponds to Vj/U Vrr/U Viv, while is the intersection of the 
rotation axis of with the d.o.c. The event horizons of the black ring and 
the spherical component coincide with ({p = 0} \ Vd. . c ) H V a . We would like 
to show that > on Vj U V77 U V/jj U V}y. Moreover, since G y vanishes 
as p 2 on the axis of [TJ, then it is necessary to check that 

lim ^ > 

p^o+ p 2 

on 

It turns out to be convenient to view as a function of p, p±, . . . , p 5 , 
where pi = \J p 2 + (z — a^ 2 — (z — a^. In this parametrization, the transla- 
tional symmetry of z, a, is explicit. The analogues of the set£| (12. ip are now 
in R 6 



V^ = {p < pi < p 5 < p 4 < P3 < P2] , 

Vj = {p>o}nv„, 

v n = {p = o,/i!>o}nv;, 

An = {p = o,/ii = o,p 5 > 0} n v,, , 

v IV = {p = o,p 4 = o,/i 3 >o}n^, 

Vp = {p = o,/i 2 = o}nv;, 



(2.2) 



and, of course, Vd. . c . C ViUt^jUVmUViyUV^. It follows from the definition 
of pi that if pi = pj for i ^ j then p = pi = pj = 0. This parametrization 
turns out to be more helpful in completing the proof. 



4 These sets are not equivalent to (|2.1|) because V4 C Viv, V2 C Vp, where Va, Vi 
correspond to the sets that are defined in the old parametrization as {p = 0, z = 04} n V a , 
{p = 0, z = C12} CiV a , respectively. 
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The numerator and the denominator of g^, when written as polynomials 
in p, pi, and some Ci, c 2 , q contain tens of thousands monomials [UI2]- One 
may check with a direct Mathematica calculation that some non-trivial 
factors from the numerator and the denominator cancel and the original 
form of may be simplified to 

_ p 4 p 5 A 2 - p 3 B 2 
9 ** " ^ ' (2 ' 3) 

where 

A = + c 2 qp 1 pzPz) + c 1 p 3 p i (-qp 2 p 5 + c 2 p 1 p i p 6 )) , 

B = (PsPi(<?PiP2 - c 2 Pip 2 p 3 ) + c 1 p 2 p 4 p i (p 1 p 5 p 5 + c 2 qp 3 P6)) , 
Pi = (>3 - Pi)(PiP4 + P 2 ) , 

P2 = (/X2 - ^4) (/il/i2 + P 2 ) (/i2/i3 + P 2 ) , / 2 4 n 

P3 = (P2 - Pl)(P2P5 + P 2 ) , 
P4 = Pl(P5 - Pi) , 

P5 = (P2 - Pl)(P2 - P4)(P2P3 + P 2 ) , 
P6 = (PlP2 + P 2 )(P2P5 + P 2 ) , 

and Pi > 0. The functions H x , F are non- negative and they were defined 
in [2]. It follows from the analysis in [1] that zeros of H X F exist only for 
p = and they cancel with the zeros of the numerator of g^p. 

The parameters ci, c 2 , g depend only on Oj and do not depend on p, z. 
However, if one assumes that p, p i are independent variables then ci, c 2 , g 
are finite continuous functions of p, pi and are given by 



(P3 - Pl)(P4 - Pl)p5(PlP3 + P 2 )(PlP4 + P 2 



PlP3P4(P5 - Pl)(PlP5 + P 



2^ 



(2.5) 



CiC 2 p 4 (p2 - Pl)(piP2 + P 2 ) 

ClPl(p 2 - P 4 )(P2P4 + P 2 ) + C 2 p 2 (p 4 - Pl)(piP4 + P 2 ) ' 

where we imposed (I2.5P on ci in the formula for g. The formula for c 2 in 
terms of p, pi is to long to be usefully cited here. It may be derived from 
equations (4.2), (5.1) in p], but it is not necessary for our calculations. 

The formulas (12. 5p . (12. 6 j) are explicitly valid in Vi U Vu- They are also 
valid in the remaining part of the d.o.c. provided the limit p — > + is carefully 
taken. If pi(p = 0) = then taking this limit should be preceded by the 
substitution pi — > p 2 fii, where pi > p]. 
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The simplification ( 12. 3 j) is a significant one. Even if one does not substi- 
tute formulas for C\, c 2 , q the original expression for written as a rational 
function contained 106995 monomials. This number was reduced to 2344 in 
the simplified formula. 

In order to present the proof in a concise form, we introduce the following 
operators. Let N, D denote operators acting on a rational functions that 
return polynomials: a numerator or a denominator, respectively. The result 
is not unique and the action of N, D is given only up to an overall factor. 
However, this non-uniqueness is not important for our problem. We also 
define substitution operators S q , respectively S c 2, that return the rational 
function which is obtained after q, respectively cf, has been substituted in 
the original expression using (12.61) . respectively (12.51) . One should note that 
S q and S c 2 do not commute with N and D in general. The formula for 
c 2 , in contrast to the formula for c±, contains a square root that cannot be 
eliminated by taking c|. We prefer to preserve the polynomial form of the 
evaluated expressions, hence we will substitute only c\ and avoid substituting 
c 2 - 

The calculations described below were done with Mathematica. We 
present them here in a brief form. It follows from the smoothness of the 
Black Saturns [I] that the expressions evaluated below are regular. 

We start with the analysis of on Vj (p > 0)@. The denominator of 
is given by p\p±H x F and it is a positive function on Vj. The numerator 
of g^ is equal to where 



If > then g^ > on V/, as desired. Firstly, we check at a random 

point@ P E Vi that E±\p > 0. Since both E + and S_ are continuous in 
p, pi, then it is sufficient to show that they cannot vanish. The functions 
S ± are linear in q, as may be seen from (12 .4p . (12 .7p . We substitute q into 
H± and examine the numerators of the resulting expressions. We would like 
to show that none of them (NS q S±) has zeros. By inspection, we find that 
NS q S± are quadratic in ci- Since C2 is real, a negative discriminant of NS q 'E± 
with respect to c 2 would imply that none of the equations NS q E± = has 
a solution. We calculate these discriminants A± and they turn out to be 

5 Vi = {p > 0} n V M as indicated in (|2.2j) . but for the sake of brevity we will remind 
only first part of the definitions. 

6 We impose the equation satisfied by C2 in such checks, but this is actually not necessary. 




(2.7) 
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fourth order in c\. Next, we substitute c\ into A± using (12. 5 p and taking 
c\ = (c 2 ) 2 , c\ = C\c\. With a help of Mathematica we have derived 

S c i&± = w±x 

/il^/i 4< U5 /2 (/il - /i 2 ) 2 (/i3 - Atl)(At 4 - /Xi)(/X 2 ~ P4) 2 (PlP 2 + P 2 ) 2 (^l/i4 + P 2 ) 2 

(A*5 -Pi)(PiP5 + P 2 ) 2 

The factors multiplying iu± are strictly positive, and iw± are complicated 
polynomials in p, y/JTi. These polynomials are linear in c\ (with non- vanishing 
coefficients in front of c\ as it will follow from our further analysis). We check 
at a random point P' G Vj- that w±\p> < 0, hence if w± have no zeros then 
A± < 0. Let cf be solutions to the equations w± = 0. A Mathematica 
calculation reveals that cf = — cj~. We set U = (cf) 2 — S c 2c\ and calculate 

4/i 2 /i 3 M/ii - A i 5) 2 (A i iA i 5 + P 2 ) 2 U ' 

where U, U are complicated polynomial^ in p, p« with signs unknown so 
far. The coefficient in front of U/U is strictly positive. Now, we succeeded 
in making the signs of U, U explicit by writing them in terms of the new 
positive functions 

A 5 i = p 5 - pi , A 45 = p 4 - p 5 , A 34 = p 3 - p 4 , A 23 = p 2 - P3 • 

The coefficients in U, U in front of p, pi, A^-, turn out@ to be positive and 
belong to 

[9, 13705432344] n Z , [1, 137075730] n Z , 

respectively. Since p, pi, A^ are strictly positive, then it follows that U, 
U are strictly positive and the equation U — does not have solutions. 
Therefore, cf ^ S c 2Cx and the polynomials w± cannot vanish. This means 
that the discriminants of NS g E± in respect to c 2 are negative (A± < 0) and 
there are no real c 2 that would satisfy any of the equations H± = 0. Finally, 
this implies that > on Vj (for p > 0), as desired. To complete the proof 
it remains to repeat the analysis above in the remaining part of the d.o.c. 

7 The polynomial U is a full square. 

8 We have /ii > p, but do not have to form monotonically increasing sequence like 
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For p = the denominator of (given by pip^HxF) is not strictly 
positive any more. However, it follows from the Black Saturns' smoothness [T] 
that whenever the denominator of vanishes the numerator of g^ (equal 
to E+EL) vanishes as well and the limit lim p ^ + is finite. Moreover, the 
continuity of g^ implies that this limit is non-negative, possibly zero. 

On Vjj (p = 0, pi > 0) the argument proceeds along the same lines as 
for Vj. Intermediate expressions have different form, but the reasoning is 
analogous. We have found 

u = $ Z. 

4At3A*4M5(A*i - A^) 2 U 

The coefficients in the polynomials U, U in front of pi, Ay range in 

[9, 7882] n Z , [9, 184] n Z , 

respectively. Hence, U > and none of the expressions NS q 'E.± vanishes, and 
g^ > on V n . 

In order to study g^ on Vm (p — pi = 0, p 5 > 0) we substitute pi = p 2 pi 
into E± (the numerator of g^ is given by S + S_). It turns out that p 4 
factors out in each term E±. On the other hand, p 8 factors in fiip^HxF (the 
denominator of g^). We set 

E± = lim —- 

P ^o+ p 4 

and repeat the proof for E± as in the case p > 0. We have 



4/i?P3Ai4/i5(l + A1P5) 2 U ' 

Finally, the coefficients in the polynomials U, U in front of the strictly posi- 
tive functions pi, A^ are again greater than zero and in 

[9, 8714] nZ, [9, 184] n Z , 

respectively. Therefore, g^ > on Vm, as expected. 

We continue our analysis on Vjy (p = Pa = 0, p 3 > 0). Here, pi = p 2 pi, 
p 5 = p 2 p5, P4 = p 2 p4- We note that < pi < p 5 < p 4 . The calculations 
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are similar to the calculations for Vjjj. A factor p 16 appears both in the 
numerator and the denominator of g^p, thus we define 

H± = lim ^| , 

p^o+ p 8 

and apply our standard analysis to E±. The formula for U is 
U = ^ U . 

4/}i/V*i/wi - h) 2 ( l + fa^) 2 

In the final step we introduce A51 = p 2 A 5 i, A 45 = p 2 A 45 and the coefficients 
in U in front of strictly positive functions pi, A 5 i, A 45 , A 34 , A 23 are greater 
than zero and belong to [9, 5852] H Z. Hence, > on V/y, as desired. 

In the case of the rotation axis of the periodic Killing field (the set 
Vrp where p = pi = 0), the analysis is slightly more involved. We substitute 
Pi = p 2 /ii into E±. Then, we have verified that E± ~ p 15 for generic ci and 
g. However, a Mathematica calculation reveals that for our choice of c 2 
and q (the equations (12. 5p . (I2.6P ) the leading terms vanish and we have at 
least S-t ~ p 16 . Therefore, we drop the leading terms in H± and analyse the 
remaining higher order terms. We denote them with 5' ± . Next, we set 

H± = lim ^| , 

p->o+ p ib 

and apply our standard procedure to E±. We have found that 

4p!p^p 4 pf (p 4 - faf(fa - faf(fa - Al)(A2 - faf 



S C? A = 



- Ml 



which is strictly negative because < p 4 < fa < p 4 < P3 < fa. Thus, we 
have ~ p 32 . The denominator of g^ (given by pip 4 if x F) behaves like 

p 30 and vanishes like p 2 . This, together with positivity of g^ for p > 
and continuity of g^, implies that 

lim 9 J!± > 

on the axis of the rotation of and the inequality (II. ip holds. It completes 
the proof of stable causality of the Black Saturns' d.o.c. 
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The generalized Weyl coordinates cover also the sets 



v br = { P = o, fj, 5 = o, /i 4 > 0} n , 
v sb = {p = 0, fi 3 = 0, p 2 > 0} n , 

which have been ignored in our analysis so far. The set Vbr = V br U V4 
corresponds to the event horizon of a black ring and the set Vsb = V sb U V 2 
correspond to the event horizon of a spherical black hole@ The closure of 
Vd.o.c. is given by Vd.o.c. U Vbh U Vsb- Since V4 C V/y, V 2 C K/,, it remains to 
apply our analysis on VJ, r , V^&. 

On Vb r (p — 115 — 0, p 4 > 0) we substitute pi = p 2 pi, ps = p 2 p5 into 
S ± . Then, it turns out that p 7 factors out of E± and p 14 factors out of the 
denominator of (p 1 p i H x F). We define 

H± = lim =^ , 

p^o+ p' 

and calculate the discriminants A± of NS q E± in respect to C2. Next, we 
substitute c 2 into A± and obtain 

Pip|pfpfp 4 (l +pi/i 2 ) 2 (l +Pl/i 4 ) 2 (/i2 -/W4) 2 
^ C 2^± = -W± " , 

where the coefficient behind w± is strictly greater than zero. This time w± 
are polynomials in p, pi, p 5 , p 4 , p 3 , p 2 . They become explicitly positive 
if written in terms of pi, A 51 , A 45 , A 34 , A 23 with coefficients in the range 
[3, 16] H Z. Then, S C »A± < and > on V br . 

The analysis of positivity of on V sb (p = p 3 = 0, p 2 > 0) mimics 
the calculations on V br . There are the following changes. We substitute 
pi = p 2 pi, P5 = p 2 p$, Pi = p 2 p^, P3 = p 3 p4 into E±. We set 



and obtain 



5 C? A ± = 

4pip 4 pj^4(l + pip 2 ) 2 (l + P3P2KI + P4P2)(1 + P5P2)(P3 ~ Pl) 3 (P4 ~ Pi) 

^5 - £l 



3 The sets V4, V2 were defined in the footnote on page 3. 
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which is strictly negative. Therefore, > on V s b- 

In summary, our analysis applied to Vbr, V s b establishes that > 
there. The strict positivity of g^ remains valid on V4 C Viv, V% C V^, as we 
know from the analysis of the sets Viv, V^. Therefore, g^ > holds on both 
event horizons Vbr, Vsb- This implies that the inequality (11. ip is satisfied 
on the event horizons. They may be included to the Black Saturns' domain 
of stable causality and the Black Saturns are stably causal on the closure of 
the d.o.c. 

We stress that we have not assigned a particular value to C2 in our calcu- 
lations so all results hold also for the Black Saturns with conical singularities 
on the axes of the rotation of the periodic Killing fields d^, d v . 
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